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1. Introduction 

The European Twisted Mass Collaboration (ETMC) has recently started an intensive, system- 
atic program of calculations of three-point correlation functions using the large number of gauge 
configurations produced for three values of the lattice spacing and various lattice volumes adopting 
the tree-level improved Symanzik gauge action and the twisted mass fermionic action with Nf = 2 
dynamical flavors tuned at maximal twist (see [|l]]). The aim is the determination of the electro- 
magnetic (e.m.) and weak semileptonic form factors relevant for light and heavy-light mesons as 
well as for baryons. 

In this contribution we present the preliminary results obtained so far for the charge form factor 
of the pion, the universal Isgur-Wise (IW) function, the vector and scalar form factors relevant for 
decays and the D ^ K{n) transition. 

The presently completed runs correspond to three simulated sea-quark masses, anisea = 0.0040, 
0.0064 and 0.0100, at jS = 3.9 corresponding to a = 0.087(1) fm (a^^ ~ 2.3 GeV) |@], and at a 
single lattice volume {V *T = 24^ *48). The mass of the spectator valence quark is fixed at the sea- 
quark mass, while the values of the mass of the valence quark struck by the electroweak current are 
taken from the set {0.0040,0.0064,0.0085,0.0100,0.0150, 0.022,0.027,0.032, 0.25,0.32,0.39, 
0.46}. The first five masses correspond to the "light" sector and coincide with the values of the 
sea-quark mass adopted by the ETMC at j3 = 3.9, the subsequent three are around the "strange" 
quark mass and the heaviest four he in the range from the "charm" quark mass nic to twice mc. At 
each value of the sea-quark mass we have computed the two- and three-point correlation functions 
for charged pseudo-scalar mesons, using the standard local 75 interpolating fields, on a set of 240 
independent gauge configurations, separated by 20 consecutive HMC trajectories. 

In order to improve the statistical accuracy we have calculated the two- and three-point corre- 
lation functions employing all-to-all quark propagators estimated through the "one-end" stochastic 
method of Ref. [|3|]. The advantages of such a procedure with respect to the "standard" one based 
on the use of point-to-all quark propagators with fixed point sources are clearly illustrated in Fig. |l]. 

Figure 1: Charge form factor of the pion, 
Fn{q^), versus the squared four-momentum 
transfer in lattice units, calculated at the 
(bare) quark mass am — 0.0040 using a subset 
of 80 ETMC gauge configurations. The dot- 
ted line is the behavior expected from vector- 
meson dominance at the simulated quark 
mass. The errors are purely statistical obtained 
by the jackknife procedure. 

Moreover, in order to get rid of the limitations in the minimum value of the spatial momentum 
imposed by periodic boundary conditions, we make use of twisted boundary conditions [Q] on 
the valence quark fields ^ We have adopted the Breit frame, where initial and final mesons have 
opposite spatial momenta, because in such a frame for a given value of the squared four-momentum 

' The use of different boundary conditions for sea and valence quarks produce finite size effects that are exponentially 
small [|]. 
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transfer q the spatial momentum injected to the active quarks is minimized. Thus one has 

2 



(1.1) 



where M, (Mf) is the initial (final) meson mass and is a real variable. In our simulations we have 
chosen d = (6,6,6) with 6 ranging from 0.1 to 1. 

Thanks to the tuning at maximal twist many physical observables at zero momentum are au- 
tomatically 0(a)-improved As for the matrix elements of the electroweak current at non- 
vanishing momenta, the improvement can be realized by a suitable averaging of matrix elements 
with meson momenta of equal magnitude but opposite sign ||^. In the Breit frame this is equiva- 
lent to the exchange of initial and final quark masses^. As a byproduct, the matrix elements of the 
e.m. current are automatically 0(a) -improved at any momenta in the Breit frame. 

We employ on the lattice the local vector current which needs to be renormalized. The renor- 
malization constant Zy can be calculated using the matrix element of the time component of the 
(local) vector current between pions at rest. Indeed, since the charge form factor of the pion is 
normalized to unity at q-^ = 0, one has 



Zv{7i+{d)\hfu-^dfd\7t+{d)) = 2M^. 



(1.2) 



Another way to obtain Zy is the use of the axial Ward Identity as carried out in Ref. The two 
determinations exhibit a quite remarkable statistical precision (~ 0.03%), and they agree very well 
in the chiral limit, as shown in Fig. ^ while at non-vanishing quark masses they differ mainly by 
terms of the order of a^mAgcD- 




Figure 2: Renormalization constant of the lo- 



cal vector current Zy calculated via Eq. (1.2) 
(full dots) and through the axial Ward iden- 
tity (squares) ^ for three values of the (bare) 
quark mass am. The errors are purely statisti- 
cal obtained by the jackknife procedure. 



2. Charge form factor of the pion 

The charge form factor of the pion, Fji{q^), is directly related to the matrix element of the time 
component of the renormalized (local) e.m. current by 

F,{q^) = ^{n+{d)\lufu-\dfd\7l+{-d)) . (2.1) 

^Note that the added correlator is almost costless if the multisolver algorithm is adopted for the inversion of the 
Dirac equation. 
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where E^ = JmI + {2kG /LY and = -A{lnd /Ly. Note that the value of is independent of 



the simulated pion mass. The matrix element appearing in Eq. (gjj) can be extracted from a suitable 
ratio of three-point to two-point correlation functions. The quality of the plateaux is illustrated in 
Fig. |[ We remind that in the two-point correlator the identification of the pion ground state starts 
already at a time around t /a= 10 (see [P]). 
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Figure 3: Ratio of three-point to two-point correlation functions, R-ipts/iptsiO^ versus the lattice time for 
two values of the (bare) quark mass am — 0.0040 (a) and 0.0060 (b), corresponding to physical pion mass 



around 300 and 380 MeV, respectively. The plateaux of Rjipts/iptsit) provide directly the form factor (2.1 ) 



The results obtained for Fji{Q^ = —q^) in the unitary setup (i.e., equal valence and sea quark 
masses) are reported in Fig. ^(a) and compared with experimental data from Ref. [10]. It can 
be seen that: i) the lattice results exhibit a remarkable statistical precision; ii) thanks to the use 
of twisted boundary conditions the form factor is precisely determined at values of as low as 
0.05 GeV^; and iii) the lattice results, obtained at pion masses of about 300,380 and 470 MeV, 
overestimate the experimental data in the whole range of values of Q^. 
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Figure 4: (a) Charge form factor of the pion Fj^iQ^) versus = —q^ (in physical units). The open markers 
are the experimental data of Ref. [p^, while the full markers are the ETMC results. The lattice points at 
Mjj ~ 380 and 470 MeV are slightly shifted in for a better reading, (b) Squared pion charge radius (in 



physical units) versus the squared pion mass (in lattice units). The experimental point is from PDG [|1 
The dotted and dashed lines are the results of ChPT fits at one-loop and two-loops, respectively (see text). 

The ^^-dependence of our ETMC results can be very nicely fitted assuming a pole behavior. 
The corresponding values of the pion charge radius are shown in Fig. ^(b) and lie well below the 



experimental value from PDG [11]. 

The predictions of Chiral Perturbation Theory (ChPT) for the pion charge radius have been 
elaborated at one-loop in Ref. [12] and at two-loop in Ref. [ p^ (in the continuum and infinite 
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volume limits). At one-loop one has 

where F is the pion decay constant in the chiral limit and Ae a low-energy constant (LEC), while 



the two-loop formula of Ref. [ ]13| ] can be rewritten as 



= HiAl/Ml) - 1] +AMl + BMl log{Ml), (2.3) 

where A and B depend on various LEC's. The value of the decay constant F has been determined 
for our unitary setup in Ref. [Q] (aF = 0.0534(6) corresponding to F ~ 121 MeV). Thus the ex- 
perimental value of the pion charge radius fix the value of the LEC Ag in the one-loop formula 
(^^), namely £^ = [log(Ag/M^)]M^=i39.6 MeV = 14.4(3), as well as its mass dependence, as shown 
in Fig. ^b) by the dotted line. The differences of the one-loop ChPT prediction and the ETMC 
points represent a clear indication of important contributions from higher loops. 

In this preliminary analysis we neglect both finite size and discretization effects, which are 



nevertheless expected to be small, and we use the two-loop formula (2.3), containing three free 
parameters, to fit both the three ETMC points and the PDG value, obtaining £(, = 17.2(7), as 
shown by the dashed line in Fig. |^b)). It can be seen that the chiral enhancement expected at low 
pion masses is hardly visible at the simulated pion masses. More lattice points, particularly below 

~ 300 MeV, are necessary to find out a clear signature of the chiral logs. 

We have calculated also the scalar form factor of the pion, limiting ourselves only to the 
connected insertion of the scalar density operator. The results will be presented elsewhere ||T^. We 
simply want to point out that the values obtained for the scalar radius exhibit the same qualitative 
features of those discussed above for the pion charge radius, including the relevance of two-loop 
effects in their pion mass dependence. 



3. Universal Isgur-Wise function 

The investigation of heavy-to-heavy e.m. transitions, described by a single form factor 



lEps 



{ps{e)\hfh\ps{-e)) , (3.1) 



allows to determine the IW function (co) by performing the infinite heavy-quark limit, viz. t,{(o) = 
lim„,,^^oo Fps{q^), where ft) = 1 — /2M^^ . Note that the form factor Fps{q^) is automatically 
normalized to unity at = for any simulated mass because of the conservation of the e.m. current. 

We have calculated Fps{q^) for various values of the (bare) heavy-quark mass, arrih, taking 
the spectator-quark mass, anisp, to be equal to the sea-quark mass arusea- In Fig. ||(a) we have 
reported the results obtained at the lowest sea-quark mass, (xnif^p — Qifi*;^^ — 0.0040. It can be seen 
that the dependence upon the heavy-quark mass is very mild so that the extrapolation to the infinite 
heavy-quark limit can be safely neglected. 



■'with such a definition we disregard the (small) perturbative correction that should be removed from Eq. (3.1) to 
get the proper renormalization-group invariant definition of (o)). 
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Figure 5: (a) The IW function calculated at am^p = anisea — 0.0040 for various values of the heavy-quark 
mass am-h. The various lines are fits based on a "pole" ansatz. (b) Values of the IW slope p/jy obtained 
through "pole" (full dots) and BSW (open squares) ansatz. The vertical line correspond to the value of the 
(bare) light u{d) quark mass determined in lattice units in Ref. 

Our results for ^ (ft)) can be nicely fitted using eitlier tlie pole or the BSW JT^ ] ansatz. Tiie 
corresponding values of the IW slope p/^y = — [d^{(o)/d(o]Q}=i are reported in Fig. ||(b). A naive 
linear extrapolation in the spectator quark mass to the (bare) light u{d) quark mass, as determined 
by ETMC in Ref. [|], gives = 0.77 ± 0.28, where the error is statistical only. To our knowledge 
this is the first determination of the IW slope with Nf = 2. Recently a quite precise value of the IW 
slope, p/^y = 0.89 ± 0. 17, has been obtained ain the quenched approximation in Ref. [p^. 

4. Kij, decays and heavy-to-light transitions 

As it is known, the matrix element of the vector weak current between pseudo-scalar mesons 
involve two form factors, the vector f+{q^) and the scalar fo{q^) ones, namely: 



{pi+P2)^f+{q ) + {pi-P2)tif-{ 

^2 



Ml-Ml 



f-{ 



(4.1) 



where = {p\ — . \n this contribution we limit ourselves to illustrate in Fig. ^ the nice level 
of statistical precision achieved in the determination of the vector and scalar form factors relevant 
for the cases of Kb decays and of the D ^ K transition. 
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Figure 6: Vector /+ {q^) and scalar /o {q^) form factors of pseudo-scalar mesons relevant for decays (a) 
and for the D ^ K transition (b). 
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5. Conclusions 

We have presented preliminary results for the electroweak form factors of light and heavy-light 
pseudo-scalar mesons, obtained at a single lattice spacing (a ~ 0.09 fm) and at a single lattice vol- 
ume (y*r = 24^*48), using the tree-level improved Symanzik gauge action and the twisted mass 
Wilson action with Nf = 2 dynamical flavors tuned at maximal twist. The use of all-to-all quark 
propagators computed with a stochastic method, as well as of twisted boundary conditions on the 
valence quark fields has allowed us to achieve both a quite remarkable statical precision and very 
low values of the four-momentum transfer. We have analyzed the mass dependence of the pion 
charge radius using Chiral Perturbation Theory, obtaining evidence of relevant two-loop contribu- 
tions. Results at more values of the sea quark mass, as well as the investigation of volume effects 
and continuum extrapolation, are however required in order to draw definite quantitative conclu- 
sions. The universal Isgur-Wise function has been computed from heavy-to-heavy electromagnetic 
transitions and its slope in the case of u{d) spectator quarks has been found to be = 0.77 ± 0.28, 
where the error is statistical only. 
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